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ABSTRACT 

The Young Stellar Object (YSO) W33A is one of the best known examples of a massive 
star still in the process of forming. Here we present Gemini North ALTAIR/NIFS laser- 
guide star adaptive-optics assisted A"-band integral-field spectroscopy of W33A and 
its inner reflection nebula. In our data we make the first detections of a rotationally- 
flattened outer envelope and fast bi-polar jet of a massive YSO at near-infrared wave- 
lengths. The predominant spectral features observed are Br 7, H2, and a combination 
of emission and absorption from CO gas. We perform a 3-D spectro-astrometric anal- 
ysis of the line emission, the first study of its kind. We find that the object's Br 7 
emission reveals evidence for a fast bi-polar jet on sub-milliarcsecond scales, which is 
aligned with the larger-scale outflow. The hybrid CO features can be explained as a 
combination of hot CO emission arising in a disk close to the central star, while cold 
CO absorption originates in the cooler outer envelope. Kinematic analysis of these 
features reveals that both structures are rotating, and consistent with being aligned 
perpendicularly to both the ionised jet and the large-scale outflow. Assuming Kep- 
lerian rotation, we find that the circumstellar disk orbits a central mass of ;>10M Q , 
while the outer envelope encloses a mass of ~15M Q . Our results suggest a scenario 
of a central star accreting material from a circumstellar disk at the centre of a cool 
extended rotating torus, while driving a fast bi-polar wind. These results therefore 
provide strong supporting evidence for the hypothesis that the formation mechanism 
for high-mass stars is qualitatively similar to that of low-mass stars. 
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1 INTRODUCTION 

It is still unclear how massive stars (M±^> 8Mq) are formed. 
This is due in part to the rapid formation timescales of such 
stars (< 10 5 yrs), which mean that the star can reach the 
main-sequence (MS) while it is still heavily embedded in its 
natal molecular cloud. Also, it is still unknown how a mas- 
sive protostar can continue to accrete despite the emmense 
outward radiation pressure it exerts on its su rroundings (for 
reviews of massive sta r formation see e.g. iBeuther et al.l 
l2007l : lHoare et alj|2007l l 

It has been suggested that massive stars may form 
through a mechanism similar to that of lower mass stars, 
whereby matter is accreted from a circumstellar disk, with 
the outward radiative force escaping via the poles and driv- 
ing a fast bipolar wind. Several numerical studies have 
succ eeded in creating stars with masses >30Mp)th is way 
fe.g. lYorke fc Sonnhaltei1l2002i ; lKmmholz et ai1l2009l '). while 



it has been shown that the fast bipolar winds should 
strongly influence t he infalling material from the envelope 
l|Parkin et all l2009h . The observational evidence to sup- 
port this scenario however is limited and circumstantial. 
The CO bandhead emission at 2.3/ mi has been modelled 
as ar i sing in a circumste l lar disk (e.g . Chandler et al.|[l993l . 
1 19951 : iBik fc Thil l2004al ; iBlum et all |2004 ): large-scale bi- 
polar outflows have bee n observed to origina te from heavily 
embedded objects (e.g. IBeuther et all 12002! ); and some au- 
thors have found ev idence of accretion disks in individua l 
massive protostars (|Patel et all 120051 ; iBeltran et all 120061 : 
iTorrelles et~alll2007l l 

Massive Young Stellar Objects (MYSOs) are thought 
to represent the phase in a massive star's formation imme- 
diately before the star reaches the main-sequence (MS) but 
whilst it is still accreting, and as such are important objects 
with which to study the formation of massive stars. W33A 
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Figure 1. Wide-field ii'-band image of W33A, taken from the UKIDSS Galactic Plane Survey |Lucas et al . 2008). The image is logarith- 
mically scaled in units of sigma above the background. The NIFS field-of-view is indicated by the white box. Our sky field is indicated by 
the black box in the lower-right of the image. The natural guide star, used in the tip-tilt correction stage of the adaptive optics system, 
is indicated by the black circle. 



is a well-known MYSO. It lies at a distance of ^3.8kpc and 
has a source luminosity of ~10 5 L m dFaundez et alll2004h . It 
is deeply embedded (e.g. iGibb et al.ll2000fi . but is visible in 
the K band despite this. The morphology of the object's 
large-scale nebula suggests a bipolar outflow, but as yet no 
convincing evidence exists for an accretion disk or a bipolar 
wind close to the central star. The object's weak radio emis- 
sion has only been marg i nally resolved at radio wavele ngths 
(|Rengaraian fc Holll996l ; Ivan der Tak fe Mentenll2005T ). The 
close environm ent of W33A has b een probed using mid-IR 
interferometry (Ide Wit et alj|2007t de Wit et al, 2009, sub- 
mitted) , with the visibilities well fitted by a model in which 
the emission arises in the warm dusty envelope close to the 
surface of the cavity created by a bipolar outflow. 

To provide a detailed study of W33A's inner nebula, 
and ultimately to investigate the physical processes occur- 
ing in the formation of the central star, in this paper we 
present high spatial resolution integral field spectroscopy of 
W33A. These data allow us to study the spatial variations in 
emission from the ionised outflow, the dense, hot molecular 



material close to the star, and the cooler molecular material 
in the object's inner envelope and outflow. 

We begin in Sect. [5] with a description of the observa- 
tions and data reduction steps. We present the results of 
the data in Sect. O and provide a discussion in Sect. [4] We 
conclude in Sect. \5\ 



2 OBSERVATIONS & DATA REDUCTION 
2.1 Observations 

In Fig. [1] we show the wide-field /('-band image of the 
W33 A region taken fro m the UKIDSS Galactic Plane Sur- 
vey (|Lucas et al.l 1200*81 ) . The image clearly shows what ap- 
pears to be a large-scale outflow extending ~0.8' towards 
the south-east at a position-angle (PA) measured to be 
(135±5)°. Diffuse emission is also seen to the north-west 
of the flux peak, which may be the opposite lobe of the out- 
flow. A velocity resolved 12 CO map of this region conclu- 
sively shows that the material is outflowing from the central 
source (see Fig. 2 of de Wit et al., submitted). 
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Our data were taken using the Gemini Near-infrared 
Integral Field Spectrograph (NIFS, iMcGregor et alj|2003t ) 
on the nights of 16 April and 25 May 2008. The observations 
were pointed such that the bright flux centre was at the 
northern edge of the instrument field-of-view, in order to 
obtain data on the base of the outflow as well as the central 
source. The instrument applies the image-slicing technique 
to separate the 3"x3" field-of-view into 29 slices, which are 
then aligned on the detector and passed through a diffraction 
grating. We employed the K grating in combination with 
the HK filter and central wavelength set at 2.2/im. This 
setup achieved a spectral resolution of R « 5300 over the 
wavelength range 2.0-2.4/im across the full NIFS field-of- 
view. 

Our observations made use of the Gemini ALTAIR 
laser-guide star (LGS) adaptive optics system. The system 
uses a natural guide star for tip-tilt correction, and the laser 
guide star for wavefront sensing. For tip-tilt correction we 
used the star 2UCAC-25155527, which has muc = 14.6 and 
is located 6.9" from the centre of our field-of-view (see Fig. 
0. 

The observational strategy involved observing the tar- 
get and neighbouring sky in an ABBA pattern, the location 
of the sky field is shown in Fig. [1] The integration time per 
individual observation was 180s, with each integration read- 
out in the faint_mode setting of 16 non-destructive reads. 
A total of 7 separate integrations were used, giving a total 
integration time of ~20mins. In addition to the science tar- 
get, we observed the AO V double HR6798 to characterize 
the telluric absorption. The GCAL calibration unit was used 
to take continuum lamp observations plus associated dark 
frames for flat-fielding purposes. Argon-Xenon arc exposures 
were obtained for wavelength calibration, as well as ronchi 
flats which are used to characterize the spatial distortion of 
NIFS data. 



2.2 Data reduction 

The preliminary stages of data-reduction were done using 
the Gemini instrument-specific pipeline which runs on the 
I RAF platform. Science observations had their associated sky 
frames subtracted to remove the dark current and sky emis- 
sion lines. The data were then divided through by the nor- 
malized flat-field to correct for pixel-to-pixel variations in 
sensitivity, and any bad pixels flagged in the earlier reduc- 
tion steps were interpolated over. The arc and ronchi flat 
frames were used to derive the transformation matrix of the 
2-dimensional 'sliced' image onto a 3-dimensional (x,y,X) 
datacube. The application of this transformation rectifies 
the data both spatially and spectrally, and it is at this step 
that the wavelength calibration is applied. 

Removal of telluric absorption and combination of the 
individual datacubes was done in IDL, using custom-written 
routines and the GSFC IDL librarjQ. Correction for the Br 7 
absorption in the telluric standard was done by dividing 
through by a synthetic AO V template spectrum, interac- 
tively altering the centroid and absorption strength of the 
Br 7 feature in the template until a satisfactory match to 
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Figure 2. Image of the AO v double HR6798A+B, illustrating the 
performance of the ALTAIR-LGS adaptive optics system during 
the observations. Contours show the encircled energy at 68% and 
95% of the total flux. The FWHM was measured as 0.13". 



the standard star was found. At each spatial pixel of the sci- 
ence datacube, the spectrum was extracted and divided by 
the telluric standard. If the spectrum had sufficient signal- 
to- noise (i.e. over ~100 counts in the continuum) the two 
spectra were first cross-correlated to achieve optimal telluric 
correction, particularly important in the region of the OH 
absorption lines at 2.0-2.1/im where sub-pixel misalignments 
can introduce artifacts into the spectrum. 

The individual science datacubes of repeat observations 
were spatially aligned by first taking the median image of 
each cube across all spectral pixels, then cross-correlating 
these median images with one another. We then took the 
sum of the spectra at each spatial pixel, using the median of 
these spectra as a template from which to reject any further 
bad pixels and cosmic-ray hits. 



2.3 Adaptive optics performance 

In order to illustrate the spatial resolution of the data and 
the PSF correction that was achieved through the ALTAIR 
adaptive optics system, in Fig. [2] we show the wavelength- 
integrated image of the telluric standard stars HR6798A+B. 
To better demonstrate the shape and size of the corrected 
PSF, the image was first split into two, with each sub-image 
containing one star. The sub-images were normalized such 
that the corrected flux from each star was the same, then the 
sub-images were recombined. On the image we have drawn 
contours indicating 67% and 95% encircled energy. While 
the star to the north appears slightly elliptical, the southern 
star is circular to within the errors. For each star we find 67% 
of the PSF flux is contained within a radius of ~0.15" of the 
flux centre. The full-width half-maximum (FWHM) of each 
star was measured to be 0.13". 
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Figure 3. NIFS K-band image of W33A when integrated over 
all spectral channels, with North up and East to the left. The 
image is scaled in units of sigma above the background, estimated 
from the NW corner of the field. The location of maximum flux 
is labelled 'A', and the other four bright knots of emission are 
labelled 'B'-'E'. 



2.4 Variations in spectral resolution and 
wavelength calibration 

The accuracy of the wavelength calibration and the consis- 
tency of spectral resolution was determined from analysis 
of datacubes constructed from the arc lamp observations. 
For each arc line that was used to determine the wavelength 
solution we fitted gaussian line profiles, and measured the 
line- width and velocity offset. We found that the absolute 
wavelength calibration had a root-mean-square standard de- 
viation of ±4 km s _1 across the field-of-view, with no dis- 
cernable trends with spatial position. The spectral resolution 
was similarly found to be Av = 55±4 km s _1 across the field, 
again having no trend with spatial position, other than that 
it was systematically poorer by ~6km s _1 in the south-east 
0.5" x0.5"region. As very little data was used from this re- 
gion, we did not correct for this effect. 



3 RESULTS 

3.1 Nebular morphology 

In Fig. |3]we show the total two-dimensional image when the 
datacube is summed across all wavelength channels. The 
central source is not a point-source, and instead is notice- 
ably extended in the north-south direction (see Fig. [3). In 
particular, we see a bright central peak (labelled 'A' in Fig. 
[3]) and a 'finger' pointing southward (labelled as 'B'). The 
'L'-shaped inner nebula, which can be marginally resolved 
in the UKIDSS if-band image and which forms the base of 
the outflow, is resolved in our data into three bright knots of 
emission, which we label 'C, 'D' and 'E'. As we show later, 
the emission from these knots is predominantly scattered 
light. 
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Figure 4. Spatial profile of the flux centroid of W33A when 
compared to the average profile of the two standard stars. While 
the profiles in the x-direction are similar, the ^/-profile of W33A 
is noticably extended in comparison to the standard. 



Table 1. Equivalent widths of the atomic emission lines. The 
uncertainty on EW was measured from fluctuations in the con- 
tinuum. 



Line 




EW (A) 


Br 7 


2.1662 


-5.00±0.17 


Mgn 


2.1374 


-0.43±0.17 


Mgn 


2.1438 


-0.24±0.17 


Nal 


2.2062 


-0.50±0.17 


Nal 


2.2090 


-0.37±0.17 



3.2 Spatial variations of spectral properties 

The unique feature of IFU data is that it is possible to chart 
the variations of the observed spectral properties across the 
field of view. We begin in this section with a description 
of the principle spectral features, followed by an analysis of 
how these features vary with spatial position. 

In Fig. [5] we plot the spectrum of the bright central 
source extracted from region A (see Fig. |3J). The steeply- 
rising continuum indicates that the central source is heavily 
reddened. In order to highlight the spectral features, in the 
bottom panel of the figure we show the spectrum when nor- 
malized by a 4th-degree polynomial fit to the continuum 
(labelled as 'region A' in the Figure). In this normalised 
spectrum we see emission features of Br 7, H2, Nal and the 
CO 2-0, 3-1, 4-2 and 5-3 bandheads. Also seen are narrow 
absorption lines in the region of the CO bandheads, at- 
tributable to low J transitions of the CO v=2-0 state, as 
well as a broad absorption feature around 2.28^im which we 
attribut e to CH3OH ice (a lso observed in this object's spec- 
trum bv lTaban et"al]l2003l ) . We do not confirm the presence 
of He 1 emission detected bv lTaban et~aH (|2003T l. while there 
is a suggestion of emission from the Mg 11 doublet at 2.14/^m. 
Tabled] shows the measured equivalent widths of the atomic 
emission lines as measured from the region A aperture. 
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Figure 5. Spectrum of region A, and the ratio of region A's spectrum to that of the average background spectrum. The top panel 
shows the spectra before continuum fitting, to illustrate the high degree of reddening. The bottom panel shows the spectra after dividing 
through by fits to the continua, to highlight the discrete spectral features. 



3.2.1 Spectrum of the central source versus surroun 
nebula 

To analyse spectral differences between region A and the rest 
of the field, we have created an average 'background' spec- 
trum by integrating over all spatial pixels in the observed 
field except those in region A. In the top panel of Fig. we 
have also plotted the ratio of the region A spectrum to this 
background spectrum, labelled 'region A/background'. The 
continuum of this ratio spectrum is shallower than that of 
the central source, which is to be expected if the emission 
in the southern nebula is predominantly scattered light. 

In the bottom panel of Fig. [5] we plot the ratio of the 
spectra when normalized by a continuum fit. Here the emis- 
sion due to Br 7, Nal, Mgll and CO disappears, suggest- 
ing that the emission morphologies of these features follow 
that of the continuum, i.e. they are also scattered. The H2 
lines go into 'absorption' indicating that this emission orig- 
inates mainly in the surrounding nebula rather than on the 
bright source at region A. The CH3OH and CO absorption 
are still seen in the residual spectrum, implying that there 
is an increased level of opacity due to these features along 
the line-of-sight to region A compared with the surrounding 
nebula. 
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Figure 6. Three-colour image of the continuum emission. Each 
RGB channel is scaled linearly, while the dynamic ranges were 
chosen arbitrarily to best illustrate the contrast between region 
A and the southern nebula. 



3.2.2 Variations in continuum emission 

In the previous section we showed that the slope of the K- 
band continuum emission is not uniform across the field. To 



further illustrate the spatial variations in continuum emis- 
sion we created narrow band continuum images at three 
featureless regions of the spectrum: 2.083/im, 2.185^im, and 
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Figure 7. Left: Image of the continuum-subtracted line flux of Br 7. The images are scaled in units of sigma above the background, 
measured from the north-west corner of the image. The overplotted contours show the continuum emission and are drawn at 20a, 30cr, 
40<T, 50<T, 70c, lOOcr and 200<r above the background. Right: Same as the left-hand figure but for the first CO bandhead at 2.293^tm. 



2.273/im, each with a bandwidth of ±0.003/im. These images 
were then combined into an RGB image, with the short- 
est wavelength image assigned to the blue channel and the 
longest to the red. The resulting image is shown in Fig. [6] 

In the image, the southern nebula is seen predominantly 
in 'blue' light (i.e. the shortest wavelength), while the central 
source appears red. This behaviour can be understood as a 
combination of scattering and extinction. The central source 
at region A is heavily absorbed, making the direct light from 
region A appear red. Meanwhile, the emission to the south 
is mostly light from region A which has been scattered, and 
so appears blue. 

3.2.3 The Bry and CO emission 

In this section we map the spatial morphology of the emis- 
sion due to Br 7 and CO (the CO absorption is studied in 
Sect. l3.2\6)l . First we subtract the continuum at each spatial 
pixel of the datacube by fitting a 4th-degree polynomial to 
featureless regions of the continuum. We then integrate the 
emission across the spectral feature in question. For the CO 
emission, we took only the bright blue edge of the first band- 
head (2.293-2.300^im) - this feature contains the majority 
of the flux, and is uncontaminated by the low-J absorption 
lines. 

The continuum-subtracted emission images of Br 7 and 
CO are shown in Fig. [7J Both images display morphology 
similar to that of the continuum image (traced by the con- 
tours), with each of the five emission knots clearly visible. 
It is unlikely that the extended Br 7 and CO emission is 
formed 'in-situ' in the southern parts of the nebula, as the 
two spectral features require different temperatures and den- 
sities. It is far more likely that in the southern nebula, as 
with the continuum emission, we are seeing scattered light 
which originates within a few AU of the central source. 

In Fig. [8] we further explore the scattered light expla- 
nation for the Br 7 and CO emission by taking the ratio of 
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Figure 8. The ratio image of the continuum-subtracted Br 7 and 
CO emission. Contours are the same as those shown in Fig. [7J 



the two continuum-subtracted images. If the emission in the 
southern nebula attributable to these transitions is due to 
scattering by material in the outflow, we should in principle 
see very little morphological difference between the emission 
morphologies. Figure[S]shows that no features corresponding 
to the southern nebula are seen in the ratio image. Instead, 
we see a slight increase in the ratio of Br 7 to CO in the 
southern regions compared that around region A. This can 
be understood as a combination of the wavelength depen- 
dence of scattering, which boosts the Br7/CO ratio in the 
south, and extinction by dense material around the central 
source, which attenuates Br 7 with respect to CO around 
region A. 
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9. The variations of the Br 7 spectral feature across the five bright knots of emission identified in Fig. [3] The mean of the five 
is overplotted as a dotted line. The dashed line indicates W33A's systemic velocity. 
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Figure 10. The variations of the CO feature across the five bright knots of emission identified in Fig. [3] 



3.2.4 Variations in line-profile with viewing- angle 

As we have shown in previous sections, the line-emission ob- 
served in the southern nebula is most likely scattered light 
which orginates at the flux centre (i.e. region A). We can 
then use this scattered light to study the emission from the 
central source as viewed by the southern nebula, when com- 
pared to viewing the emission directly from region A. In 
effect, this allows us to study the line-profiles as a function 
of viewing angle. This is a technique we have previously 
used to study th e wind geometry o f the evolved massive 
star IRC +10420 (|Davies et al.ll2007l) . 

We explore the spatial variations of Br 7 and CO in 
Figs. [5] and 1101 where we plot the integrated spectra of the 
five bright knots of emission identified in Fig. [3] as A-E. 
In respect to the Br 7 emission, there is very little varia- 
tion across the five regions. The regions show evidence for 
P Cygni absorption, and perhaps an extended red wing. The 
equivalent widths (EWs) of the profiles in the southern neb- 
ula are slightly larger than those of regions A and B. As 
the Br 7 emission in the southern nebula is most likely scat- 
tered light, these small changes in EW are possibly due to 
anisotropy in the ionized circumstellar gas. 

The CO emission shows clear variations across the ob- 
served field. To highlight this in Fig. [11] we plot the spectra 
in the vicinity of the v=2-0 bandhead. The blue edge appears 
shallower in regions A and B, while the edge is much steeper 
in the three southern regions. Also, the transitions between 
2. 30-2. 305pm appear to be narrower and more pronounced 
in the southern knots (C, D and E) than in those close to 



region A. This suggests that the CO emission, when viewed 
directly through the bright central source at region A, has 
a broadened velocity profile with respect to the profile re- 
flected from the southern nebula. As the southern nebula's 
viewing angle to the central source is different from our di- 
rect line of sight to the same region, this suggests that the 
CO profile's broadening is a function of the viewing angle. 
Such behaviou r is typically attributed to the emiss ion form- 
ing in a disk l|Bik fc Thill2004al : iBlum et al.ll2004h . As the 
material in the southern nebula 'sees' less rotational broad- 
ening than our direct line-of-sight, the orientation of the 
disk is consistent with being aligned perpendicularly to the 
larger-scale outflow. 



3.2.5 Spectro-astrometry of Brj 

Here we investigate the small-scale geometry of the Br 7 
emission by performing a 3-D spectro-astrometric analysis 
of the Br 7 line. At each wavelength channel of the dat- 
acube we extracted a 2-D image and determined the astro- 
metric centroid of the emission at region A. This method 
can give spatial information on extremely small scales - the 
astrometric accuracy on the centroid is roughly equivalent 
to A ~ FWHM/SNR. So, for the spatial resolution of our 
observations (0.13") and a SNR at region A of several hun- 
dred, in principle we can expect to achieve accuracies of 
order of a tenth of a milliarcsecond. In order to verify this 
we performed spectro-astrometric tests on the data from 
the telluric standards. We found that indeed the r.m.s. on 
the centroid position was of order O.lmas. We tried different 
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Figure 11. Spectra of the five regions identified in Fig. \3\ around the CO v=2-0 bandhead. For clarity, spectra have been interpolated 
onto a finer grid. It is clear from this figure that the line profiles of C, D and E are narrower and more pronounced than those of A and 
B. 
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Figure 12. Spectro-astrometry across the Br 7 emission line. (Left): the bottom panel shows the intensity spectrum of region A, while 
the top panel shows the astrometric position of the centroid flux peak as a function of wavelength. The red line, labelled 'DEC, shows 
the centroid position in the y-direction, while the black line ('RA') shows the position in the x-direction. (Right): the same data as in 
the left-hand panel, but represented in x — y space. The symbols are colour-coded to show the velocity at each point. The dotted line 
shows the position-angle of the red and blue high-velocity components. 



methods of centroiding, including flux-weighted centroiding, 
finding the location at which the derivatives go to zero and 
fitting 2-D gaussian profiles. All techniques were found to 
give similar results within the errors. In addition we observed 
no artifacts around the telluric absorption lines, and no ex- 
cursions across the telluric standard star's intrinsic Br 7 ab- 
sorption, as is to be expected (see Appendix A). 

In the left panel of Fig. [12] we plot the position of 
the centroid as a function of wavelength in the region of 
Br 7. The spectro-astrometric profiles have had any large- 
scale gradients removed by fitting a 3rd-degree polynomial 
to the continuum. Firstly, we note that the astrometric ac- 



curacy of the centroid, as judged by the continuum regions, 
is comparable to that predicted (0.1-0.2mas), corresponding 
to a physical scale of 0.4-0.8AU for a systemic distance of 
3.8kpc. Secondly, the centroid shows a clear shift in the re- 
gion of Br 7. In the north-south direction (labelled 'DEC 
in the figure) , the centroid appears to shift south slightly in 
the region of the P Cygni absorption, before shifting north- 
wards by ~1 mas. Similarly, in the east- west direction (la- 
belled 'RA', where the positive direction is eastward), there 
is a slight shift to the east coincident with the P Cygni ab- 
sorption, followed by an excursion to the west of ~0.5mas. 
This behaviour is illustrated again in the right panel of 
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Fig. 1121 This time, the centroid is plotted in RA-DEC space 
with colour-coded symbols to indicate the velocity of each 
point. The location of the continuum centroid is illustrated 
by the cross. The plot shows that the emission with the 
largest red-shifted velocities (~400km s _1 ) is located to the 
north-west of the continuum location, while the emission 
blue-shifted by a similar amount is located to the south- 
east of the continuum, indicating expansion. From analysis 
of those points with speeds between 300-400 km s" 1 we find 
that the statistical significance of the red- and blue-shifted 
points being spatially separated from the continuum is 2.8a 
and 4.1ct respectively. In addition to this behaviour there 
is a large excursion across the centre of the line with \v\ < 
100 km s" 1 . 

The observed behaviour can be explained if we make 
the assumption that the Br 7 emission line is made up of 
two or more components. The broad high- velocity emission 
originates in an axisymmetric structure, most likely a bi- 
polar jet, while the narrow low velocity emission is formed 
in an extended structure such as a small, dense Hn-region. 
The nature of this narrow emission is difficult to disentan- 
gle, but it is possible that it is formed in an aspherical 
clumpy medium with much velocity-dependent extinction 
an d self-absorption. This is similar to the p icture deduced 
bv lDrew etail l|l993l ) and lBunn et al.l (|l995l ) from velocity- 
resolved line-profile ratios. 

In addition to the detection of an outflow, the spectro- 
astrometric data also indicate that the bright flux maximum 
observed in the integrated image is indeed the location of the 
central star. 

For the broad-line component of the emission we can 
measure the PA of this structure by taking the data at 
300 < ln/kms -1 ] < 400 and taking the mean positions of 
the red- and blue-shifted components. We measure the PA 
of this structure to be 113±25° (see right panel of Fig. I12[) . 
and so is aligned with the large-scale outflow seen in Fig. [T] 
(PA=135±5°). From the alignment of these angles it is likely 
that in this small-scale structure we are seeing the base of 
a bi-polar ionized wind or jet. In addition, the blue-shifted 
portion of the Br 7 emission is oriented towards the south- 
east, consistent with the hypothesis that the outflow seen 
in the wide-field image is the lobe coming towards us, while 
the red-shifted lobe to the north-west is largely obscured. 

The outflow velocity, as measured from Fig. 1121 appears 
to be around ti out j ow cosi ~ 300 km s~ . Assuming an incli- 
nation of 60° (de Wit et al., 2009 submitted), this gives an 
outflow velocity measured from Br 7 of ~600km s -1 . This 
is likely an underestimate of the outflow's terminal veloc- 
ity, as a large fraction of the Br 7 emission may form in 
the accelarating part of the wind. This value is consistent 
with M YSO jet velocities measured by radio proper motion 
studies (|Marti et al.lll998l ; ICuriel et al.ll2006l ). 

3.2.6 The CO absorption 

The low-J transistions of the CO v=2-0 state are observed 
in absorption in W33A's spectrum (Sect. 15^)) . The popula- 
tion of these states requires much cooler conditions (~30K) 
than that of the high- J states seen in emission (fe w x 1000 K) 
l|Chandler et all 1 19931 . Il995l ; iBik fc Thil l2004al ). Therefore 
these two distinct spectral features must arise in different 
regions - the CO bandhead emission likely originates in ma- 
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Figure 13. Image of the datacube integrated over the CO ab- 
sorption lines. The image is scaled in units of sigma below the 
continuum flux level, where sigma is measured from the southern 
nebula. 



terial close to the central star; while the absorption lines are 
formed in the outer envelope at several thousand AU (see 
later) . 

The right-hand panel of Fig. [9] shows that the CO ab- 
sorption lines vary in strength between the bright central 
source (A) and the southern regions (C, D and E). Specifi- 
cally, the absorption is stronger towards region A than the 
southern nebula. To explore this behaviour further, sub- 
cubes were extracted from the primary data-cube, each con- 
taining the full spatial information of one of these absorption 
lines. We observed a behaviour which was similar for each 
line and so all lines were combined to increase the signal-to- 
noise. 

In Fig. [13] we show the image obtained when the dat- 
acube is integrated over these features. The image shows 
that the CO absorption is concentrated around region A 
in a band extending roughly east- west, and is perpendicu- 
lar to the large-scale outflow. Significantly less absorption is 
seen in the southern nebula, indicating that the material in 
the southern regions sees less obscuring material than our 
direct line-of-sight to region A. This suggests that the enve- 
lope has the morphology of a torus, with a physical radius of 
~2000 AU, and with an opening to the south through which 
the outflow escapes. 

Figure Q3] shows the measured velocity centroid of the 
CO absorption at each spatial position across the extinction 
structure seen in Fig. 1131 The image has been smoothed by 
a boxcar filter of width 5 pixels to highlight any velocity 
gradient. At negative (blue-shifted) velocities, the opacity 
appears to be concentrated in a region about 0.2-0.3" east 
of the centre of the bright source (i.e. region A). At positive 
(red-shifted) velocities however, there is a clear shift in the 
location of the opacity centre to a region 0.1-0.2" west of 
region A. No apparent trends with velocity are seen in the 
southern regions of the nebula, though the data are very 
noisy and are not shown in the figure. We determined the 
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Figure 14. The velocity centroid of the CO absorption lines as 
a function of spatial position across the field. Velocities are indi- 
cated only for those absorption lines which have strengths greater 
than 5cr below the continuum. The data has been smoothed with a 
boxcar filter of width 5 pixels (0.25") to accentuate the regions of 
maximum and minimum velocity. Dotted lines indicate the rough 
position angle of the structure. 

PA of rotation to be 60±10° by fitting the blue and red 
flux peaks in Fig. 1141 This is therefore consistent with being 
perpedicular to the larger scale outflow of PA=135±5°. The 
PA of this feature is perpendicular to the large scale outflow 
observed in Fig. [1] as well as to the small-scale outflow seen 
in the spectro-astrometric analysis of the Br 7 emission. As 
with the broadening of the CO bandhead emission described 
in the previous section, the observed behaviour is indicative 
of an extended rotating structure seen edge-on. The kine- 
matics of this feature are explored further in Sect. 3] 

3.2.7 Molecular hydrogen 

The light of H2 2.12/^m shows a very different morphology 
to that of the continuum and the Br 7 and CO emission 
(see Fig. I15|) . The emission is concentrated to the south-east 
of region A, while very little H2 emission is seen at region 
A itself. The other H2 lines observed in the spectrum dis- 
play similar behaviour, though have much poorer signal-to- 
noise and so are not shown. The orientation of this emission, 
extending roughly south-eastwards of region A, is approx- 
imately coincident with that of the base of the large-scale 
outflow seen in the wide-field image of Fig. [T] 



4 DISCUSSION - CO KINEMATICS, AND THE 
MASS OF THE CENTRAL OBJECT 

The pieces of evidence presented above are all consistent 
with the qualitative picture of a massive forming star in 
which the central object is surrounded by a rotating disk 
from which it is accreting. A rotationally-flattened envelope 
surrounds the central star and its disk, while the star drives 
a fast bi-polar wind. 

In this section we present a kinematic analysis of the 
CO bandhead emission and the low-J CO absorption lines 



H 2 flux 




Offset (arcsec) 

Figure 15. Same as Fig. [7] but for H2 2.12/im. The point-like 
emission at the very centre of region A is possibly due to poor 
continuum subtraction. 



and use them to derive the mass of the object about which 
they rotate. 



4.1 The bandhead broadening 

As shown in Fig. 1111 the CO bandhead emission lines appear 
to be broader towards region 'A' than those observed in the 
southern nebula. This behaviour can be attributed to a flat- 
tened rotating structure (i.e. a disk) which is seen at lower 
inclination when viewed directly than when seen reflected by 
the nebula to the south. If we assume the extreme case that 
(a) the line-of-sight 'A' views the disk edge-on, and (b) that 
the southern nebula views the disk pole-on, by convolving 
the profile seen from the southern nebula until it matches 
that at 'A' we can estimate a lower limit to the disk rotation 
speed. By assuming Keplerian rotation we can then estimate 
the mass of the central object about which the disk rotates. 

In Fig. [16] we plot the CO bandhead profile seen at 'A'. 
We also plot the profile seen in the southern nebula (region 
'C') which has been convolved with a rotational broaden- 
ing function of the form ip(v) ~ [1 — (v/Kot) 2 ] 1 ^ 2 , where 
v is velocity and V ro t is the disk rotational velocity. The 
best match to the blue-edge of the bandhead, as well as 
the lower-J line profiles, is found when the rotational veloc- 
ity is 50±5km s _1 . If the direct line-of-sight is not exactly 
edge-on, and/or if the orientation of the reflected light is not 
exactly pole-on, the disk rotation speed will be greater than 
this value. By taking the inclination angle recently derived 
by de Wit et al. (submitted) of 60±20° (i.e. 30° away from 
being edge-on), we can further constrain the disk rotation 
speed to 641 2 , 2 km s _1 . 

Assuming that the disk rotation is Keplerian, and that 
the CO emission arises in a narrow annulus at a fixed radius 
R from the central object we can now derive a mass of the 
central object M c from, 

M c = RV 2 /G (1) 
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Figure 16. The observed profile of the CO bandhead in region 'A' (black solid line), as well as the profile seen in the southern nebula 
(from region 'C') when convolved with a PSF of various widths. 



where G is the gravitational constant. iBik fc Thl l|2004bj ) 
who found that the CO typically forms at a distance from 
the central star of 1- 3AU. Using R = 2±1 AU, and a distance 
to W33A of 3.8kpc (|Faundez et alj|2004h . we find that the 
lower limit to the central mass is 



4.2 The absorption-line velocity structure 

In Sect. l3.2T6l we showed that the absorption lines of CO ap- 
pear to be formed in a rotating structure coincident with re- 
gion A, and measured the PA of this structure to be 60±10°. 
In Fig. [17] we show a position- velocity slice across the aver- 
aged CO datacube at a PA of 60° . The data shows the clas- 
sical signature of a Keplerian rotating structure. To help 
illustrate this, overplotted are models of Keplerian rotation 
for three differe nt central masses and with an assumed dis- 
tance of 3.8kpc (|Faundez et alj|2004h . Each curve has been 
convolved with the spatial resolution of our observations 
(0.13"). A formal least-squares fit to the data yields a cen- 
tral mass of M sin i = 13±2A/0 , a result which is insensitive 
to the uncertainty in PA. If we again assume the system in- 
clination is 60±20° we find that the outer envelope sees a 
central mass of 1S^_\Mq. Note that the central mass seen 
by the envelope will include the mass of the envelope itself. 
This number is very similar to that derived for the central 
mass enclosed by the circumstellar disk, and suggests that 
the majority of the system's mass is contained within the 
central star. 

If the central object has arrived on the main-sequence, 
this mass would correspond to a spectral type earlier than 
Bl and a ze ro-age min-sequence lum inosity of log(L/ Lq = 
4.3±0.3 fe.g. lMevnet fc Maederl200o| ). This is in good agree- 
ment with the latest estimation of W33A's bolometric lumi- 
nosity from MIPS photometry, log(L/L0= 4.7 (Mottram 
et al. 2009, in prep). The derived spectral-type and lumi- 
nosity may also explain the detection of the Mgn emission 
lines. Such lines are typically seen in the spectra of early- 
to-mid B supergiant s, and are primarily formed in the wind 
jHanson et al.l l2005). The ionization emission lines that we 
see in W33A (Br 7 and Mg n) are probably formed in the bi- 
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Figure 17. Position-velocity diagram across the rotating struc- 
ture seen in the CO absorption lines. The image is scaled in terms 
of fractional intensity below the continuum. The crosses indicate 
the velocity of the intensity minimum for the absorption as a 
function of positional offset, and contours are drawn at intensity 
levels of 0.895 to 0.92 at intervals of 0.005. Lines show models 
of Keplerian rotation for three different central masses and an 
assumed distance of 3.8kpc. 

polar outflow, the density of which may be similar to that 
of a typical B supergiant wind. 



5 SUMMARY AND CONCLUSIONS 

We have presented a study of adaptive-optics assisted near- 
infrared integral-field spectroscopy of the young stellar ob- 
ject W33A. The results can be summarized as follows: 

• In the first 3-D spectro-astrometric study of its kind, we 
have resolved the Br 7 emission on sub-milliarcsecond scales 
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(physical scales <1 AU), and find that the emission geome- 
try is that of a fast bipolar wind which is aligned with the 
large-scale outflow seen in wide-field images. The Hb 2.12/jm 
emission also traces an outflow with the same orientation out 
to 

• The extended nebulosity to the south is seen in the 
continuum, as well as in the light of Br 7 and CO. We have 
shown that this emission is light reflected by the outflow 
lobe, and so can be used as a 'mirror' with which to observe 
the emission from the central object from different angles. 

• We interpret the CO emission as arising in a disk which 
has a low inclination with respect to our line-of-sight. In the 
southern nebula we see light reflected towards us from the 
disk at a lower inclination angle (i.e. more pole-on), which 
causes the bandhead profile to be less velocity-broadened. 
This is consistent with the disk plane being perpendicular 
to the large-scale outflow. The velocity broadnening implies 
a lower-limit to the central mass of M c = 10tjjM Q . 

• From the low-J v=2-0 transitions of CO at ~2.35/im 
we find evidence for a rotationally-flattened cool molecular 
envelope, or 'torus', at a radius of ~2000 AU from the central 
star. The plane of rotation is perpendicular to the small- and 
large-scale outflows. Analysis of the velocity structure of the 
torus indicates that it is rotating about the outflow's axis, 
and that it is orbiting a central mass of 15^M Q . 

Our findings suggest a picture of massive star forma- 
tion within W33A which is in excellent qualitative agree- 
ment with the accretion-disk-plus-bipolar wind paradigm. 
An accretion disk orbits a massive central star, which is 
surrounded by a cool molecular envelope which has been 
rotationally-flattened into a torus. The central star is driv- 
ing a bipolar wind, seen on small scales in the ionized gas, 
and on larger scales as molecular / continuum emission. The 
two measurements of the central mass indicate that at most 
the accretion disk makes up ~30% of the system, and so the 
system mass is dominated by the central star. 
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Figure Al. Spectro-astrometry of the Br 7 line for the telluric 
standard star used in our observations, similar to Fig. 1 12 1 The 
lack of artifacts across the telluric features and the star's Br 7 
absorption, as well as the observed level of astrometric precision 
(~0.2mas), validate the results of the similar analysis of W33A. 

Yorke, H. W. & Sonnhalter, C. 2002, ApJ, 569, 846 



APPENDIX A: SPECTRO-ASTROMETRY OF 
THE TELLURIC STANDARD STAR 

In Fig. lAll we show a spectro-astrometric analysis of the tel- 
luric standard star following a similar methodology to that 
presented in Sect. 13.2.51 The intensity spectrum (bottom 
panel) shows the Br 7 absorption line, as well as the many 
telluric features. The top panel shows the astrometry of the 
star's flux peak as a function of wavelength, illustrating that 
no discernable variations are seen across neither the Br 7 
line nor the telluric lines. This of course is to be expected, 
and the Figure serves to show that the technique we employ 
produces no spurious artifacts greater than the noise level 
(~0.2mas) which could be misinterpreted. 



